Symmetric and non-symmetric vortex-antivortex molecules 
in fourfold superconducting geometry 
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In submicron superconducting squares in a homogeneous magnetic field, Ginzburg-Landau theory 
may admit solutions of the vortex-antivortex type, conforming with the symmetry of the sample 
[Chibotaru et al, Nature 408, 833 (2000)]. Here we show that these fascinating, but never experi- 
mentally observed states, can be enforced by artificial fourfold pinning, with their diagnostic features 
enhanced by orders of magnitude. The second-order nucleation of vortex-antivortex molecules can 
be driven either by temperature or applied magnetic field, with stable asymmetric vortex-antivortex 
equilibria found on its path. 

PACS numbers: 74.20.De, 45.05+x, 74.78.Na, 74.25. Dw 



All sub-atomic particles have a fraternal twin: an an- 
tiparticle, exactly alike except for e.g. opposite charge or 
helicity. However, as a general rule, matter and antimat- 
ter cannot coexist without excess energy and annihilate 
each other. 

This universal duality has its analogue in the physics of 
superconductors, where vortices as the carriers of mag- 
netic flux play the role of charged particles. Namely, 
vortex-antivortex pairs in superconducting films can be 
induced in a local hot-spot, created by thermal fluctu- 
ations [H, photon absorption or driving current 
Intuitively, vortex-antivortex pairs are easily stabilized 
in an inhomogeneous magnetic field, such as one result- 
ing from dipolar magnetic objects in the vicinity of the 
superconductor 4]. 

Contrary to latter examples, recent theoretical stud- 
ies P, 6, 7] have shown that vortex-antivortex structures 
can be stabilized in submicron superconductors even in 
homogeneous magnetic field, and without any apparent 
energy input. These findings are restricted to the cases 
when vortex structure does not conform with the sam- 
ple geometry. A typical example is the C3 symmetry of 
the three-vortex state in a superconducting square; in 
spite of the unipolar applied field, the C4 state of four 
vortices with a central antivortex may become energeti- 
cally preferable. Similarly to the solid-fluid and the fer- 
romagnetic symmetry-breaking processes, the more sym- 
metrical phase is found on the high-temperature side 
of the transition, and the less symmetrical one on the 
low-temperature side. Landau pointed out that one can 
always unequivocally determine whether or not given 
state possesses a given symmetry Therefore, it is 
not possible to analytically deform a state in one phase 
into a phase possessing truly different symmetry. This 
means, for example, that it is impossible for the solid- 
liquid phase boundary to end in a critical point like the 
liquid-gas boundary. Nevertheless, the above described 
vortex-antivortex nucleation during the C3 to C4 tran- 
sition in mesoscopic superconductors is of second-order 
with respect to temperature. The appearance and sta- 



bility of these asymmetric vortex-antivortex molecules in 
the ground state is one of the main objectives of this 
Letter. 

It should be noted here that the symmetry-induced 
vortex-antivortex configurations remain experimentally 
undetected, mainly because of their high sensitivity to 
defects in sample edges ^ . As another drawback, due to 
extreme vortex proximity those states are undistin- 
guishablc from a single multiquanta vortex for conven- 
tional techniques such as scanning-tunnelling and Hall- 
probe microscopy. Theory suggests that the latter prob- 
lem can be solved by a magnetic dot deposited on the 
sample [3] . Added inhomogeneous magnetic potential in- 
teracts dually with present vortex and antivortex [lo| and 
separates them in a controllable fashion, but disturbs the 
conceptual novelty of no external stabilizing factor added 
to the vortex-antivortex coexistence. 

To address above issues, we discuss in this Letter the 
properties and the improved observation conditions of a 
vortex-antivortex state in a square superconducting sam- 
ple with fourfold pinning. We introduce a cluster of 2 x 2 
nanoholes to enforce the C4 symmetry, the non-energetic 
source of vortex-antivortex nucleation (see Fig. [T]). Al- 
beit, induced states do not always obey the imposed sym- 
metry. We show that similar symmetry breaking can be 
achieved in a classical cluster of oppositely charged par- 
ticles in a strongly confined geometry. The latter setup 
has recently emerged as the 'standard model' for a va- 




FIG. 1: Superconducting square (size w x w and thickness d) 
with four nanoholes (each Wh x Wh in size, displaced by Ih 
from X and y symmetry lines). Indicated direction of applied 
homogeneous magnetic field H is denoted as positive. 
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riety of systems on different energy and length scales, 
not only flux lines in superconductors, but also vortices 
in Bose-Einstein condensates, electrons on liquid helium, 
colloidal suspensions, and dust particles in plasmas. 

Our theoretical approach relies upon the Ginzburg- 
Landau (GL) theory, where we solve self-consistently a 
set of mean field differential equations for the order pa- 
rameter tp and the vector potential A 
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The latter is the Maxwell- Ampere equation with a cur- 
rent density j = — \ip\'^A. k is the mate- 
rial property and equals the ratio between the mag- 
netic field penetration depth A and coherence length ^. 
Solution of Eqs. lU ^ minimizes Gibbs free energy 
g/Go = J[~\Tp\'^ + 2{A-Ao)-j]dn, where the con- 
dition j± — was used on the boundaries of the super- 
conducting volume 17, and Go stands for the supercon- 
ducting condensation energy H^/Att. In above expres- 
sions, V X Aq = H denotes the applied magnetic field, 
and V X A = h the resulting local field in the sample. 
All distances are expressed in units of = C(T=0), the 
vector potential in (/)o/27r^Oj and the order parameter in 
\J —aj (5 with a, /? being the GL coefficients. For details 
of the numerics we refer to Ref. ■ 

Symmetry-induced asymmetry Fig. [2] shows the 
ground-state phase diagram of our perforated sample, 
coinpared to the previously studied one without holes 
@, IE 0, B [m ■ Namely, at given temperature T (scaled 
to critical temperature Tc), we obtain stable solutions of 
Eqs. IlllS]) by starting the iterative procedure from ran- 
domly generated initial conditions, while applied mag- 
netic field H is swept up/down. The ground state is 
then determined by comparing the energy of all found 
states. Note that dimensions of the system are given in 
^0, so that depending on the material our results apply 
to samples from hundred nanometers to several microns 
in size. Parameter k is taken equal 1. 

The size-field phase diagram [size scaled to ^(T) ~ 
Co / ^ (^/^c)^] for a square sample without holes has 
been studied previously in Refs. 1, [3|- Since in Ref. 

the linearized GL theory is used, valid only extremely 
close to the superconducting/normal (S/N) phase bound- 
ary, those results were quantitatively corrected in Ref. 
where both GL equations were numerically solved. 
However, authors kept parameter Ate// = K^i{T)/d fixed 
when temperature was changed, so their diagram did not 
reflect true temperature dependence of the vortex state. 
Nevertheless, our phase diagram (with temperature de- 
pendence taken explicitly) shows similar qualitative be- 
havior - sets of individual vortices are found at lower tem- 
peratures, which merge into a giant vortex at higher tem- 
peratures for all vorticities (L) except L — i. The trian- 
gular symmetry of that state [Fig. Wis)] is not favorable 




FIG. 2: Temperature versus magnetic field equilibrium phase 
diagram. Solid lines illustrate 1"' order, and dashed 2"'* order 
transitions in the ground state. Parameters used: k = 1, 
d = 1^0, w = 10^0, Ih/w = Wh/w = 0.125 (see Fig. [TJ. 
Contourplots (a)-(f) show the phase of the order parameter 
of particular vortex states (as indicated in the left diagram). 



closer to Tc when square confinement dominates, and the 
L=4-l(=3) vortex-antivortex state is induced [Fig. [DJa)]. 

However, contrary to previous works, our analysis 
shows that the vortex-antivortex (VAV) nucleation pro- 
cess is of second order (dashed lines in Fig. [2]) over a 
broad temperature and field range. In other words, with 
changing applied magnetic field (A(/) « (/)o) or tempera- 
ture (AT K, O.lTc) while in the L = 3 state, besides the 
three existing vortices a vortex-antivortex pair gradually 
dissociates. A remarkable non-symmetric L — A—l state 
is induced [Fig. [H^b)], and the ground-state exhibits sym- 
metry breaking during the transition between intrinsically 
symmetric states (C3 vs. C4). 

Although these asymmetric equilibria seem fairly coun- 
terintuitive, their stability follows from the competing in- 
teractions in this system. In analogy, we recall the prop- 
erties of small confined clusters of particles, where the 
confinement would result from the screening currents in 
our system (along the sample edges), and charged parti- 
cles represent vortices. 

Analogy to classical systems The model sys- 
tem was defined in Ref. where Hamilto- 
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with TV being the number of particles. We took 
square-parabolic confinement potential of size Wc, i-e. 
i„ 



V{x,y) 



1 



5 -\- cos^ (2 arctan ^ , 



(m- 



mass of particle, wg-confinement frequency) with a 
nonzero i5 ^ 1, ensuring the existence of the deriva- 
tives of this potential in the corners of the square. To 
further translate this system to vortices in mesoscopic 
superconductors, one has to choose properly the vortex- 
(anti)vortex interaction energy (Uij). Firstly, this inter- 
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FIG. 3: Energy of the 3-particle cluster compared to the clus- 
ter of 4 particles and one antiparticle, as a function of the size 
of the square confinement Wc, for modified logarithmic inter- 
particle interaction (of strength /3, and with a = 1). (a-c) 
are snapshots of the ground state configuration of the latter 
cluster with decreasing Wc {x,y in units of 



2/3 



action must be of long range (and diminishing far from 
the source), knowing that supercurrents decay ^ l/r 
away from the vortex. Secondly, Uij may not diverge 
for rij < a but saturate, allowing for realistic merging of 
vortices into a giant vortex [l^, or a vortex-antivortex 
annihilation. In this scenario, a roughly corresponds to 
the finite size of the vortex core. 

Though any interaction energy of described profile 
would suit our analysis, we used a modified logarith- 
mic and modified Coulomb interaction, more specifi- 
cally dU/dr{rij) = (3{\- e^-'^^l^- f jr^ and dU/dr{nj) = 
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)^/r^j, respectively. Notice that the mod- 



ified Coulomb interaction is a quality fit for = 
(up to a multiplying constant) to the realistic vortex- 
vortex interaction in type-II samples (see e.g. Ref. [l6j). 

Next we employed molecular dynamics simulation, a 
relatively simple and rapidly convergent technique that 
gives a reliable estimate of the energy for small clusters 
of interacting particles. The obtained energy versus con- 
finement size (wc) is shown in Fig. [31 for a cluster con- 
taining three particles of the same charge (i.e. Uij > 0) 
and a cluster of four particles and an antiparticle, with 
clear analogy to our L = 3 and L = 4 — 1 vortex state. 
For large Wc, the ground states of the two clusters have 
identical energy, as one particle and the antiparticle sit on 
top of each other. However, for tighter confinement, the 
N — 4+l~ cluster attains lower energy than the N = 3'^ 
one. Namely, particle and antiparticle gradually sepa- 
rate causing the rearrangement of the remaining 3 par- 
ticles. Figs. [Sja-c) show snapshots of this dynamical 
transition in the ground state. To our knowledge, this 
symmetry breaking in the crystallization of confined sys- 
tems is novel and applies to classical particles interacting 
with charged impurities in the substrate [IT] , and poten- 



tially to small ionic crystals under pressure [ig]. Note 
also that colloidal molecules have recently been realized 
experimentally on hydrophilic square templates [l^. 

Distinction of antivortex in fourfold pinning 

As shown in Fig. [21 in our sample with nanoholes, the 
symmetrically placed perforations (i) pin (attract) vor- 
tices individually, i.e. favor the multivortex state, and 
(ii) stabilize vortex states commensurate with the square 



geometry, e.g. L = 2 and particularly L = 4 [13|. Al- 
though L ~ 3 state is somewhat suppressed in the H — T 
space compared to the plain square case, the antivortex 
in the L=4-l state is clearly pronounced [see Fig. [lUd,e)] 
since remaining 4 vortices are captured by the holes. 
To illustrate better this issue, we calculated the vortex- 
antivortex distance {d^av) and amplitudes of the Cooper- 
pair density and local magnetic field between them, as a 
function of the sample parameters. As shown in Fig. [H 
we are able to achieve the vortex-antivortex distance of 
5-20% of the sample size, which greatly exceeds the dis- 
tances found in plain squares (see Ref. Q), especially in 
larger samples. At the same time, the amplitudes of both 
order parameter and local magnetic field are enhanced 
by almost two orders of magnitude compared to previous 
studies, and reach the limits of conventional scanning 
techniques. Our study shows that = l^ = 12.5%'w are 




FIG. 4: The properties of the L — 4—1 vortex-antivortex state 
in the perforated sample, at T = 0.84Tc, and applied fiux = 
5.5<^o: (a) the distance, maximal l^*]^ -density, and difi'erence 
in the magnetic response between vortex and antivortex, as 
a function of the position of the holes; (b) the Cooper-pair 
density and the magnetic field profile (Ah = h — H) along the 
diagonal of the sample, in case when vortices are not residing 
in holes [lu = Q.lQhw, see (a)]. 



4 



optimal parameters for the observation of the antivortex, 
as larger holes capture more applied flux and decrease the 
field in the center of the sample, whereas larger distance 
between the holes results in a weaker influence on the 
vortex-antivortex state in the sample center (see Fig. U). 
It should also be emphasized that fourfold distribution 
of holes is more important than the actual shape of the 
sample, i.e. even a circular disk with four symmetrically 
arranged holes exhibits the L = A — 1 state. 

One should note that the inter-vortex interaction is 
governed by k, and may even change sign in type-I sam- 
ples This strongly affects the vortex configura- 
tions, and consequently the vortex-antivortex stabiliza- 
tion. In our calculation for both squares with and with- 
out nanoholes, decrease of k did not favor the L = A ~ 1 
state but instead the L = 3 giant vortex becomes more 
stable, contrary to the findings of Ref. [lH. Apparently, 
suggested vortex-antivortex repulsion at lower k is dom- 
inated by attraction between the remaining vortices. On 
the other hand, for larger k, the shorter range of inter- 
vortex interactions reduces the stability of asymmetric 
states which rely on strong competing forces. Therefore, 
material properties are very important for potential ex- 
perimental observation of the antivortex, as k can be 
arguably controlled by impurities, which increase A and 
decrease ^. 

Enhanced stamina of vortex-antivortex Yet an- 
other crucial experimental issue is the sensitiveness of 
the vortex-antivortex state to surface defects. While the 
ground-state strongly depends on symmetry, the analysis 
of Ref. Q has shown that defects as small as O.Olw at the 
edges of a plain superconducting square disable vortex- 
antivortex nucleation. We found that 2x2 nanoholes 
significantly strengthen the L — 4—1 state. Fig. [5]shows 
the free energy diagrams of (i) a plain square with a small 
edge defect, (ii) a square with nanoholes and an edge de- 
fect, and (iii) square with nanoholes and a defect at one 
of the holes. Whereas in case (i) no vortex-antivortex 
state is found, in both cases (ii) and (iii) wc found sta- 
ble L = A — 1 state. Our numerical experiment shows 
that vortex-antivortex states in latter cases can survive 
for defects in edges up to 10% of the sample size and 
defects in nanoholes up to remarkable 40% of their size. 
Note also that the presence of defects changes the current 
profile of the sample, which again results in asymmetric 
vortex-antivortex configurations [see Fig. [5ljb,c)]. 

Therefore, we may conclude that artificial fourfold pin- 
ning vigorously enhances (in all relevant aspects) the ex- 
perimental observability of vortex-antivortex molecules 
in flat superconducting samples. The verification of the 
found asymmetric vortex-antivortex equilibria is of par- 
ticular importance, as this symmetry-breaking is pre- 
dicted both for many-body and few-body systems with 
strong competing interactions in polygonal constraints. 

This work was supported by the Flemish Science 
Foundation (FWO-Vl), the Belgian Science Policy, the 




FIG. 5: Gibbs free energy for T = 0.78Tc and (i) plain square 
sample with an edge-defect (size O.S^o), (ii) sample with 2x2 
holes and an edge-defect, and (iii) sample with 2x2 holes and a 
defect in hole-edges. Insets (a-c) show the Cooper-pair density 
plots of found states with total vorticity 3, in cases (i-iii). 



JSPS/ESF-NES program, and the ESF-AQDJJ network. 
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